Detection of acceleration radiation in a Bose-Einstein condensate 
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We propose and study methods for detecting the Unruh effect in a Bose-Einstein condensate. The 
Bogoliubov vacuum of a Bose-Einstein condensate is used here to simulate a scalar field-theory, and 
accelerated atom dots or optical lattices as means for detecting phonon radiation due to accelera- 
tion effects. We study Unruh's effect for linear acceleration and circular acceleration. In particular, 
we study the dispersive effects of the Bogoliubov spectrum on the ideal case of exact thermaliza- 
tion. Our results suggest that Unruh's acceleration radiation can be tested using current accessible 
experimental methods. 



One of the surprising fundamental consequences of 
relativistic quantum field theory is the dependence of 
the concept of particle number on the observer's state 
of motion. While inertial observers see the vacuum as 
empty, non-inertial observers generally perceive this vac- 
uum as populated with particles. Unruh jl| showed 
that a uniformly accelerated particle detector perceive 
the field in vacuum as a thermal state with temperature 
^■bT\j = ha/2irc, where a is the proper acceleration. The 
Unruh effect is related to other particle creation effects 
in curved space-time, such as Hawking radiation, and 
the Gibbons-Hawking thermalization in a cosmological 
expansion [2j. 

Numerous experimental ideas for detecting the effect 
have been suggested. They include, accelerated electrons 
in circular high energy accelerators circular motion of 
electrons in a Paul trapQ, intense laser induced electron 
acceleration [j| and passage of atoms through a cavity @. 
Other setups simulate the Gibbon-Hawking cosmological 
expansion thermalization effect in an expanding Bose- 
Einstein condensate 

(BEC)0, 

and in an expanding linear 

ion trap [of. (See also [Iol|). 

In this letter we propose to simulate and detect the Un- 
ruh effect using accelerated atom dots (AD) ll| or using 
optical lattices in a BEC. Since the relevant velocity is the 
speed of sound, c s « 1 [mm/ sec], Tjj « 10[nK sec 2 / m] x 
a[m/sec 2 ] and the currently feasible acceleration of op- 
tical lattices may reach a as 5 x 10 5 [m/sec 2 ], the Unruh 
temperature can be significantly higher than the relevant 
energy scales, the AD minimal energy gap (« 100Hz ss 
nK), and the BEC temperature. 

Let us begin by recalling some features of the Unruh 
effect. A detector is modeled as a localized system with 
internal levels \g) and |e) = <r + [g) and energy gap uid, 
which moves along a trajectory Xd(t) and t(r), where r 
is the detector's proper time. In the simplest free 



scalar field cj>, initially in its vacuum state, couples with 
the detector through 



Hi = g 



4>(x D (r),t(r)). (1) 



By evaluating the transition amplitudes between the lev- 
els, it is then found that for inertial trajectories the de- 
tector remains unexcited, while for uniformly accelerated 
trajectories the detector becomes thermalized. This can 
be seen by evaluating to the lowest order in g the tran- 
sition amplitudes. Inserting xd(t) — ^-cosh^f, and 
t(r) = ~ sinh 2£ , and the expression for a free field c/>(x, t) 
in Eq. ([T]), one finds that a field mode uj has a time de- 
pendent coupling of the form: g e (r,u>) = exp(i^e -i! f). 
This readily yields transition probabilities which satisfy 
Excitation/ ' Pde- excitation = e~ E l hBTv , where T v is the 
Unruh temperature. 

It is important to note that: 

i. The appearance of the effective coupling g e (r,oj) 
is sufficient in order to thermalize the detector. A 
similar coupling is also a landmark of the Hawk- 
ing and cosmological thermalization effects. 

ii. In the Unruh effect property i. is a direct con- 
sequence of the detector's accelerated motion. This can 
be easily seen 13( by noticing that the field mode uj is 
Doppler shifted in the detector's rest frame to oj'(t) = 
u)q 1 ~ v / c — ujoe~ aT / c . Therefore, the relevant col- 

lectcd phase factor becomes exp(i J w(r)rfr)) = g e (u),r). 
in. The Unruh effect is manifestly relativistic. Hence 
the interaction (PJ is defined in the detector's rest frame, 
and the trajectory, XDR.{t) = c \/t 2 + c? /a? coincides 
with non-relativistic acceleration only for sufficiently 
short times. 

The above points quantify, with increasing refinement, 
important aspects of the Unruh effect, which one wishes 
to simulate in a specific model. For example, i. can be 
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obtained by modifying the vacuum normal mode frequen- 
cies u> to uj{t) — Lue^ 2 ^, and realized in an ion traps by 
changing the trap frequency *9], or by an an expanding 
BEC [7( . In what follows we suggest a model that incor- 
porates properties i. and ii., and finally shortly discuss 
possible realizations of Hi.. 

It is well known that small perturbations of the BEC 
Schrodinger field satisfy a relativistic-like Klein-Gordon 
equation with the speed of sound c s playing the role of 
c[8(. Nevertheless, the transformation laws for a mov- 
ing detector will remain non-relativistic. We can there- 
fore obtain the effective coupling constant (i.) as a 
consequence of non-relativistic Doppler shift by choos- 
ing a modified trajectory: XDeff(t) = (c s t + ^e~ a '/ Cs ^ 
which differs from the relativistic trajectory XDFt(t) 
above (when c = c s ), by 0[a 2 t 3 /c s ] for short times, 
and 0[c 2 s /a 2 t] for long times. The Doppler shift uj' — 
luq(1 — v/c s ) = u)Qe~ at l Cs , has the same time dependence 
as in the relativistic case, with r — > t. We hence expect 
that a suitable detector that moves along XDeff will be 
similarly thermalized. 

Consider then a setup with hyperfine levels, a and d, 
where a forms a condensate described by the field \&. 
Level d will be used for an AD produced by a localized 
potential Vd ll|, [lj] or by an optical lattice. It will be 
sufficient to consider only one level with a wavefuction 
ipd{x) and creation and destruction operators d, dJ. Since 
Vd affects only atoms in the state d, in the absence of 
further coupling with the condensate, moving about Vd 
will not disturb the condensate state. We need however 
to make sure that nonadiabatic excitations are negligible. 
The adiabatic condition in this case can be derived by 
transforming to the AD rest frame and for the trajectory 
XDeffit) is given by: vT = a 2 /{ui 2 c s ) <C xq where xo 
is the width of the wave function. For trap frequency 
lu pa 100kHz, a/c s w uu d ~ 100 Hz w 10~ 3 w, and since 
the l.h.s of the inequality is less than lA, the condition 
is satisfied. Atomic levels then couple through elastic 
collisions, which to the lowest order redefin the detuning 
6, and produce self interaction terms gddd^d^dd. A large 
gdd is used 11] to simulate a two-level detector (Eq. []}. 
In the following we found more convenient to assume 
small gdd, hence the detector is a harmonic oscillator. 

We couple between the AD and the BEC by laser in- 
duced Raman transitions described by interaction Hamil- 
tonian 



are essential to the effect. However our interest is in 
the resulting coupling with phonons. Using Bogoliubov's 
theory we expand the field operator 



Hi, 



8<&d + fi, a dxifid(x) ( d^(x) + h.c 



(2) 



*(a;) = <f>(x) + u k {x)e- lult c k + v k {x)e +jMt J 



(3) 



where Vl a is the Rabi frequency. At first sight Eq. @ 
lacks the number non-conserving terms of Eq. {1} , which 



where (f>(x) is a c-number, and u k (x), v k (x) and c k are the 
phonon mode functions and annihilation operators. This 
brings the BEC Hamiltonian to a free field form Hb ec — 
J2k ^kc\c k , and spectrum uj k = (c s k) 2 + {^) 2 that is 
"relativistic", uj ~ k, for k < k c — mc s /h. 

Inserting Eq. ([3]) into Eq. ([2]), and assuming that 
ipd(x) extends over scales smaller then the phonon wave- 
length, (the dominant coupling arises from long wave- 
lengths), we obtain 

H in t = SdU+ ^VL a {d + d)) + (4) 
fid f ^2(u k (x D )e^ kt c k + v k {x D )e Wkt cl k ) + h.c, 

k 

where n a is the effective number of condensate atoms 
at the AD. For k <C k c , u k w v kl this model coincides 
with Unruh's detector model Eq. (1), apart from the 
term y^ri^n a (d + d 1i ) which describes the interaction with 
the mean-field. This term can be eliminated using a two 
mode condensate with levels a and b that couple as in 
(Eq. 2) via Raman transitions and with Rabi frequencies 
satisfying £l a = — f2&. Cancelation of this term is then 
obtained from the symmetry of the Hamiltonian. Alter- 
natively, one can use a single mode condensate and re- 
move the displacement in the AD final state by applying 
the unitary exp(-^^^(c? — d))). This approach requires 
a precise control of • x /n^[12j . 

Consider the effect of Hint on the AD when the conden- 
sate is in its ground state: c k , a \BEC) = 0. For a uniform 
motion x = vt, the excitation amplitude is to first or- 
der i £ T dtJ2 k v k {x D {t))e l ^ + ^ v ^ t ¥^ £ fe S(u d + 
(1 — v/c s )oj k ). Therefore as long as u < c s , the detec- 
tor remains unexcited. For the suggested non-inertial 
trajectory Xu e ffi as l° n g as & < v ki x ) ~ u k{x) ~ 
exp(ikx), the transition amplitudes reduce to A±(u>) oc 
J_ T exp(±iujdbt — iuj k e~ at )dt, which coincides with Un- 
ruh's expressions, with t replacing r. The total transition 
probability, is then P± = J2 k l^±(^)| 2 - F° r a given mode 
k the contribution to A±{k) comes from the saddle point 
at t ~ t s (u>) = ^loguj/uid ± T^gg i ^ or l° n g er interac- 
tion times, in order to recover the Unruh effect larger 
momenta are required, with k that grows exponentially 
with the duration of the interaction. Consequently, in a 
realistic situation the finitcness of k c will cause deviations 
for acceleration time T > t s . To resolve this difficulty we 
have to restrict the interaction time and study in more 
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FIG. 1: l(a)Thermalization of an oscillator coupled of the 
smearing time 7. The field is modeled by 20 modes, a = 
2,u>d — 1 and the coupling is f/50. The Green line denotes 
the numerical results. The Blue line is derived by numerical 
integration of the transition amplitudes. The inset shows the 
thermalization curve for 7 = 17. The red line shows the the- 
oretical temperature. The blue and green graphs shows the 
cooling and the heating curves. The results are in the units of 
ujd- 1(b) In order to check whether the distribution is thermal 
we changed the energy gap and calculated the temperature. 



detail the deviations due to the k dependence of ^ , and 
the dispersion relation. 

In a experiment of finite time T, it is important to 
consider with care the modification due to the tem- 
poral change of the coupling strength. In the sim- 
plest case of abrupt change in the coupling, the tran- 
sition probabilities can be approximated as P±(u>) ~ 
2 

A±(w) + J_ e T^dt s M _ since different modes con- 
tribute at different times, the total contribution is ef- 
fectively averaged and P± oc |^4±| 2 + -\. The correction 
does not decrease with the energy gaps since |^4±| 2 scale 
as jjL. More generally we shall assume that the coupling 
starts and ends smoothly over a time scale 7 by adding 
a regulator e _t / 7 , i.e., a slow decoupling function. 

In the following we have assumed that the detector is 
accelerated for time T and that the experiment is re- 
peated ri times by moving the AD back and forth in 
the BEC. We have studied numerically a BEC with a 
finite number of N phonon modes and described the de- 
tector by a harmonic oscillator. The total state is then 
described by a covariance matrix, and detector's popula- 
tion and temperature are derived from the AD reduced 
covariance matrix. 
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FIG. 2: The final temperature using the full Boguliubov the- 
ory, as a function of the interaction time to (with 7 = 0). The 
ratio of the energy gap and the cutoff energy was taken as 
1 : 500, hence deviation are expected for lo > 120^- This 
corresponding to acceleration times to ~ l/alogl20 ~ 2.4. 
Indeed, the theoretically predicted temperature, T ~ 0.67, is 
obtained in this simulation for t < to- 



We first considered the ideal case with k c — > oo. As is 
shown in Fig. Q^a) the effective temperature of the detec- 
tor changes gradually until it reaches a final steady state 
after n ~ 100 repetitions. The temperature is slightly 
higher then the value of Tjj since the finite decoupling 
time and the final coupling strength increas the average 
final energy of the steady state. By increasing 7 we can 
get closer to the theoretical value of T\j. In Fig. [TJb) 
the final temperature is plotted for various values of the 
detector energy gaps. As can be seen the temperature 
remains unchanged in agreement with a thermal distri- 
bution, up to fluctuations of AT/T ~ 1%. The correc- 
tions observed here is due to the finiteness of the number 
of modes and the interaction time. 

Next we extended the analysis to the full problem 
with a finite cutoff scale u) c — c s k c which corresponds 
in a realistic BEC to more than 10kHz, and is two 
orders of magnitude larger than the atomdot's mini- 
mal energy gap, which is limited by the fluctuations 
of the laser. There are two types of corrections. The 
first type is due to the changing dispersion relation; 
since the phase in the transition amplitude is now 

given by e ± iu >dt e i(kx-uit) _ e ±iui d t e i(ckt-cke~ at /a-u(k)t) _ 

e ±iuj d t e z{(ck~uj(k))t-cke at /a)^ ^ detector's energy gap is 
corrected by ck — ui(k), which is always a negative quan- 
tity. For certain modes the effective detector gap can 
vanish, which implies a divergence in the resulting par- 
tial excitation probability. For higher modes the temper- 
ature can then becomes negative, which causes a gradual 
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FIG. 3: Simulation of thermalization and finite size effects 
in circular motion. The final analytically computed tempera- 
tures is 3[fia>d] (red lines) and 5.5[hu>d] (blue lines). To exam- 
ine finite size effects we compare between two cases: the upper 
graphs for which the system size L > c a t (L/c s = 125[27r/o;d]) 
and boundary effects are minor, and the lower graphs for 
which L < c 3 t (L/c s = 25[2-7r /u>d\) and oscillations due to 
finite size effects become noticeable 

population inversion since P+ = ? 7rc ° i, m , 1 , — -, the ra- 
tio for large frequency tends to unity. This cutoff effect 
would be felt once LUd — ck — u)(k), which is smaller than 
the field cutoff. The second type of correction comes from 
the modified momentum dependence of mode functions 
Mfe and Ufe. As k increases decreases to zero, hence for 
T > t c (k c ) the temperature starts decreasing. 

In order to observe the Unruh effect, we can reduce 
the effects of the above 'ultra-high' frequency correc- 
tions by selecting a sufficiently short time scale. Fig. 
[2] displays the resulting final temperature for a numeri- 
cal computation which includes all Boguliubov's theory 
corrections. The expected thermalization effect can be 
observed but due to the shorter interaction time requires 
a slightly higher number of repetitions, n ~ 300. In or- 
der to decrease the number of repetitions the initial state 
can be chosen at the vicinity of the final temperature. To 
avoid finite temperature corrections, we need to have the 
phonon number in the relevant interacting modes to be 
smaller then 1. For a BEC temperature of 50nk this 
requires a gap energy on the order of 1kHz. 

Another interesting experimental possibility is to sim- 
ulate the effect of circular acceleration 0, [l5| . Unlike the 
ideal Unruh effect, here the accelerating detector sees 
the vacuum as excited but usually is not thermalized, 
i.e., it's final temperature depends on the energy gap. 
The advantage of this setup is that the detector does not 
have to satisfy relativistic equations and thus no special 
path is needed. Moreover, in the limit where the fre- 
quency of rotation is much smaller than the energy gap, 
the interaction is effectively non-zero only with a finite 
band of frequencies. Consequently the effect can be in- 
sensitive to the cutoff. The limit of v ~ c is especially 



interesting since the temperature divergence and the de- 
tector becomes thermalized, making this regime ideal for 
experiment. The rotation of the detector could be real- 
ized either using dipole traps [3| or optical lattices [13], in 
both setups the speed of rotation could reach the speed of 
sound. This makes the circular variant simpler to man- 
ifest than the linear one. Fig[3] displays the numeric re- 
sults for a thermalization effect of a circulating atomic 
quantum dot. 

We remark that a fuller relativistic-like realization, can 
be done as follows: we consider the condensate coupled to 
the AD as in Eq. (2), but choose the detuning (5 = 0. Us- 
ing Boguliubov's expansion we then obtain the Hamilto- 
niantf,. w n(t)y>(d+rf)+(d-rf) J2 k u k (b k -b\)] = VtH 1 ', 
where k < k c was assumed. The first term represents 
the free detector Hamiltonian (energy levels have be- 
come superpositions of number states) and the second 
term the interaction with the field. The idea is then to 
use the common factor Q(t) and modify the laser inten- 
sity so that f2(i) oc 4r = , 1 Upon integration 

J H r dt = J H'dr, hence this recovers the Unruh effect 
for a uniform-like accelerating trajectory. 

In conclusion, we found that a moving AD or an atomic 
lattice in a condensate can be used to detect acceleration 
radiation effects that are analogous to the Unruh effect. 
Our results indicate that the measurability of such effects 
is within reach of current methods. We hope that the 
analogy that we are making may be also useful the other 
way around; that is to interpret what happens when one 
moves a particle in a condensate with some acceleration. 
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